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ABSTRACT 

The conjecture that the anc i ent gl obular clusters (GCs) formed at the center of their own dark matter halos 
was first proposed by Peebles ( 1984), and has recently been revived to explain the puzzling abundance patterns 
observed within many GCs. In this paper we demonstrate that the outer stellar density profile of isolated GCs is 
very sensitive to the presence of an extended dark halo. The GCs NGC 2419, located at 90 kpc from the center 
of our Galaxy, and MGC1, located at ~ 200 kpc from the center of M31, are ideal laboratories for testing the 
scenario that GCs formed at the centers of massive dark halos. Comparing analytic models to observations of 
these GCs, we conclude that these GCs cannot be embedded within dark halos with a virial mass greater than 
10 6 Mq, or, equivalently, the dark matter halo mass-to-stellar mass ratio must be Mbm/M* < 1. If these GCs 
have indeed orbited within weak tidal fields throughout their lifetimes, then these limits imply that these GCs 
did not form within their own dark halos. Recent observations of an extended stellar halo in the GC NGC 1851 
are also interpreted in the context of our analytic models. Implications of these results for the formation of GCs 
are briefly discussed. 

Subject headings: Galaxy: globular clusters — globular clusters: general 



1. INTRODUCTION 

Despite decades of intense theoretical effort, the forma- 
tion of the ancient g lobular c lusters (GCs) remains a largely 
unsolved problem. Peebles (1984) considered the possibil- 
ity that GCs form within their own dark matter (DM) halos 
at high redshift. The growing evidence for significant self- 
enrichment in GCs and the broad acceptance of hierarchical 
structure formation has deepened interest in this formation 
scenario. Evidence against this scenario was found in the 
observations of thin t idal tails surrounding m any GCs (e.g., 
iGrillmair etaLlll995t lOdenkirchen et all 120031) . because nu- 
merical simulations showed that such t idal tails do not form 
if GCs reside within extended halos (lMoord["l996l) . How- 
ever, later work highlighted the fact that even if Milky Way 
(MW) GCs were once embedded within massive dark ha- 
los, these halos would have been tidally stripped away by the 
present epoch (Bromm & Clarke 2002; Mashchenko & Sills 
2005). This requires relatively strong tidal fields, which sug- 
gests that GCs in the outer halo of the MW may still be em- 
bedded within dark halos, if they formed within them. 

Other theories for the formation of G Cs do not appeal to 
formation at the center of dark halos. iFall & Reesl d 19851) 
proposed that GCs form from thermal instabilities in the hot 
gaseous halos expected to surround massive galaxies today. 
This proposal suffers from the fact that many galaxies that 
host GCs are not expected to reside in halos massive enough 
to support a ho t halo, such as dwarf spheroidals. 

iGunrj (T1980) was the first to suggest that GCs could form 
in the gas compressed by strong shocks. This proposal re- 
ceived tentative confirmation with the discovery of many 
massive you ng star clusters within the in t eracting Anten- 
nae s ystem (|Whitmore & Schweizeri 11995b IWhitmore et alJ 
1999) and the discovery of super sta r clusters within nearby 
galaxies (e.g., iHoltzman et al.l [T992I) . This scenario, modi- 
fied to include as formation sites any massive, dense, cold 
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patch of gas, is now the prevailing paradigm for GC for- 
mation (e.g., lHarris & Pudritzl 119941) . and, when incorpo- 
rated into our broader theory of cosmological structure 
formation, is capable of explaining a variety of observa- 
tions (e.g.. lAshman & Ze pfll 19921: iKr avtsov & Gnedinl 120051: 
iMuratov & Gnedinll20ldlh 

This prevailing paradigm for GC forma tion is complicated 
by the existence of nucle ar star clusters dBoker et al.l 120041 : 
IWalcher et al.l 12005. 2006), which implies that at least some 
GC-like systems can form at the centers of massive dark ha- 
los. The existence of young nuclear star clusters makes this 
point particularly compelling, since these clusters could not 
have migrated to the center via dynamical friction. Thus, 
while dark halos are not necessarily required for GC forma- 
tion, the conditions for GC formation may sometimes be re- 
alized at the centers of dark halos. Clearly, further constraints 
on the formation sites of GCs is desirable. 

In a series of papers, Spitzer and collaborators derived 
the kinematic properties of stars in the stellar halo of a GC , 
where stars are only ma rginally bound ( Spit zer & Hartll 19711 : 
Spitzer & Shapiro| [l972l) . An important result from this work 
was that the density profile of stars in the stellar halo should 
scale as r" 3 5 . In the present work we build upon these results 
by investigating the sensitivity of the stellar density profile to 
the presence of a massive dark halo. 

2. THE STELLAR HALOS OF GLOBULAR CLUSTERS 
2.1. Analytic model 

In this section we derive the outer stellar density profile 
of GCs embedded in a massive dark halo. The following 
derivation closely follows the assumptions and approxima- 
tions made in a series of papers by Spitzer and collaborator s 
(ISoitzer & HardlTlmllSpitzer & Shapiroll97llSpitzeill98l . 
to which the reader is referred for details. 

The density profile of a stellar system can be derived from 
its distribution function, /, via: 

«(r)oc / /(£,./) 27rv,dv,dv,., (1) 
Je<o 



2 



CONROY, LOEB & SPERGEL 



where v, and v,- are the tangential and radial velocities. We 
assume that GC halo stars are on radial orbits, and thus are 
justified in making the approximation that v 2 = v 2 , and we 
can substitute v, = J jr. Most importantly, we assume that 
f(E,J) = \E\g(J), where g is some function of angular mo- 
mentum. This functional form arises when the orbital ener- 
gies are only slightly below zero, the number of stars in the 
syste m is large, and th e system has reached a steady state 
(see lSpitzer & Shapiro! 1 19721 for details). These constraints 
require that the two-body relaxation time is short compared to 
the age of the Universe. We then have: 



«(r) oc r"Y (/) / |£|dv, 



(2) 



E<0 



where g' is some new function of angular momentum. As- 
suming that J is not a function of r in the stellar halo, we drop 
all reference to J from here on. 

For a purely stellar system we have E = iv 2 + 5>*, where $* 
is the potential of the stars and is approximated by a Keplerian 
potential ($„ oc -GM*/V). Upon substitution into Equation 
2 we recover the familiar result that n(r) oc r in the halo 
of GCs. This re sult has been confirme d by direct TV-body 
simulations (e.g., Baumgardt et al. 2002). 

Our task here is simply to re-evaluate this integral with the 
addition of a DM potential, $dm- The distribution function of 
weakly-bound stars is unchanged with the addition of a dark 
halo since the derivation makes no reference to the form of 
the potential. We therefore have: 



n(r) oc r 



£<() 
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DM 



dv, 



which upon integration becomes: 

«(r)ocr- 2 ($» + $ DM ) 3/2 . 
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We assume an NFW density profile for the dark halo that is 
motivated by collisionless ACDM cosmological simulations 
iNavarro et al.l d!996l 119971) . The implied dark halo potential 
is 

ln(l + r/r s ) 

$dm = ~G M DM g(c) , (5) 

r 

where Mum is the total dark halo 'virial' mass, c is the con- 
centration defined as c = r v /r s where r v is the virial radius and 
r s is the scale radius, and g{c) = [ln(l +c)-c/(l +c)] _1 . Over 
the physically relevant range of 2 < c < 10, g(c) varies from 
2.3 to 0.7. 

Finally then, we have the following expression for the stel- 
lar density profile in the presence of a dark halc0: 



n(r) oc r 



-3.5 



1 + 



M, 



g(c)ln(l + r/r s ) 



-,3/2 



(6) 



For Mdm/A^* <C 1 we recover the familiar result of n(r) oc 
r -3 5 When the dark halo mass is significant, the profile 
can be decomposed into three regimes. At sufficiently small 
scales the first term in brackets in equation [6] dominates over 
the second, and the profile scales as r" 3 5 . At larger scales, the 

3 The contribution from unbound stars is not included here. We expect 
their contribution to be negligible because simulations consistently find that 
stars are unbound at a rate of ~ 1 % per relaxation time. Moreover, the den- 
sity profile of the escapers is approximately r~ 2 iSpitzer 1987), even in the 
presence of a dark halo, and so their presence would not impact our conclu- 
sions. 
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FIG. 1 . — Stellar density profiles normalized to the density at 20 pc. Models 
are shown for several values of the dark halo-to-stellar mass ratio, Mdm/M* 
(top panel) and dark halo scale radius, r s (bottom panel). In the top panel 
r s = 250 pc, and in the bottom panel Mdm/M* = 10 2 . The blue and red 
dashed lines have logarithmic slopes of —3.5 and -2.0, respectively. 

second term dominates, and it takes on two limits for r smaller 
or larger than r s . For r < r s the second term scales as r and 
the total density profile then scales as n(r) oc r" 2 . At scales 
greater than r s the second term in brackets shallows, and the 
resulting density profile consequently steepens. 

In FigureQ]we show the expected stellar density profiles for 
several values of the parameters Mbm/M* and r s . For sim- 
plicity, we have fixed the virial radius to r v = 1 kpc although 
the models are insensitive to this simplification. Notice the 
strong sensitivity to Mdm/M* and the weak sensitivity to the 
r s over the scales of interest. The weak sensitivity to r s is due 
to the fact that the logarithmic slope of the dark halo potential 
varies slowly across r s . 

FigureQDdemonstrates that the density profile over the range 
10 i5 r S 100 pc is very sensitive to the presence of a dark 
halo. Our derivation of the density profile is strictly appro- 
priate only for the stellar halo of a GC, and so the profiles in 
Figure Q] will not represent real GCs on smaller scales. We 
have also ignored tidal stripping and the fact that the relax- 
ation time at large scales may under certain circumstances be 
longer than the age of the Universe. 

The derivation of the stellar density profile in the halo of 
GCs relies on the assumption that the two-body relaxation 
time is short compared to the age of the Universe. We now 
verify under what conditions this assumption is valid. The 
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FIG. 2. — Relaxation time in units of the half-mass relaxation time, as a 
function of clustocentric distance. Results are shown for several values of the 
dark matter-to-stellar mass ratio. The age of the Universe is shown in units 
of the half-mass relaxation time of MGC 1 (dotted line). 

radial dependence of the effective relax ation time in the stel- 
lar ha lo can be estimated as follows (see Lightman & Shapiro 
119781 for details). The relaxation time, t r , scales as 



t r oc 



D( AE 2 ) ' 



(7) 



where D is the diffusion rate and E is the energy. The diffu- 
sion rate is simply t he energy change per unit t ime, which, for 
stars in the halo, is (Binney & Tremaine 1987): 



D(AE ) oc 



(8) 



where P is the orbital period and e is the (small) change in 
energy per orbit. The key feature of stars in the halo is that 
they are on radial orbits that pass through the central regions 
of the GC. This fact implies that e is approximately contant 
for stars in the halo; i.e., the change in energy per orbit does 
not depend on the apocentric distance of the orbit. Thus: 



t r (xE z P. 



(9) 



Assuming that the potential is a combination of a Keplerian 
and an NFW dark halo, as we have throughout this section, 
we arrive at the following expression for the relaxation time: 



t r oc 



1 



, , M DM , , , , r/r s 

1 + — g(c) \a{\ + r/r s )- 



l + r/r s 
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(10) 



The relaxation time is a function of radius and dark matter- 
to-stellar mass ratio. These dependencies are illustrated in 
Figure [2] In this figure the relaxation time is scaled to the 
relaxation time at the half-mass radius. Notice first the coun- 
terintuitive result that in systems dominated by a Keplerian 
potential the relaxation time in the halo is actually a decreas- 
ing fu nction of radius. As discussed in Lightman & Shapiro 
(1978), this arises because the change in energy per orbit, e, 



10' r 



10- 



£ 10" 



10"- r 



irr 



- 1 1 — i — i — i — r-| 

M DM /M,= 
M DM /M» = 10° 
M DM /M» = 10' 
M DM /M, = 10 2 
M nM /M, = 10 4 



O 



• MGC1 
t- O NGC 2419 



<» *« 

o 

o • 



_l I I I 1_ 



10 



100 



Rp (pc) 



FIG. 3. — Stellar surface density profiles normalized to the surface density at 
20 pc. Our models, which include a stellar component of mass M t embedded 
within a dark halo of mass Mdm are shown as lines for a range of mass ratios. 
These models are compared to data from the GC MGC1 located in the outer 
halo of M31 (Mackey et al. 2010) and the GC NGC 2419 located in the outer 
halo of the MW gellazzini 2007D- Data are only plotted for R p > R h . The 
blue dashed line has a logarithmic slope of —2.5 and is the predicted surface 
density profile for a pure stellar system. 

is constant, while the energy of a star scales as r~ l . At greater 
clustocentric distances stars therefore require fewer orbits to 
change E 2 by of order itself. The increasing period with in- 
creasing distance is not sufficient to counteract this trend. 

The addition of a dark halo modifies this behavior, such 
that larger dark matter contributions result in longer relax- 
ation times. For sufficiently large dark matter fractions, the 
relaxation time will exceed the age of the Universe. Figure 
[2] includes an upper bound provided by the age of the Uni- 
verse in units of the half-mass relaxation time of one GC we 
will consider in the next section, MGC1. For this cluster, 
the relaxation time in the halo does not exceed the age of the 
Universe at < 6R h , for M DM /M* = 100 and at R < 30R h for 
Mdm/M* = 10. The very outer stellar halo of isolated GCs 
m ust be interpreted with these facts in mind. 

Lightman & Shapiro ( 1978) provided a simple derivation of 
the stellar density profile in the GC halo. In a steady state the 
net stellar flux through a spherical shell of radius r must be 
constant, which implies that nir)^ jt r = const. This formula 
for n(r), when combined with Equation [10] reproduces the 
stellar density profile derived earlier in this section (Equation 
|6| in the limit where E ~ 0. This is not surprising because the 
principal assumption in both derivations is the existence of a 
steady state in the halo. 

2.2. Results 

As mentioned in the Introduction, most ancient GCs are on 
orbits that would likely have resulted in severe stripping of an 
extended dark halo, were they originally embedded in such 
halos. GCs at large galactocentric distance, in contrast, orbit 
within very weak tidal fields, and so one may expect these 
objects to have retained their dark halos, if they ever had them. 

Two GCs are particularly noteworthy in this regard: NGC 
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2419 in the MW and MGC1 in M31. NGC 2419 resides at 
90 kpc from the center of our Galaxy, has a half-mass and 
tidal radius of 20 pc and 230 pc, respectively, and a V-band 
luminosity of 5 x 10 5 L© (Harris 1996b. which implies a total 
stellar mass of ss 10 6 M©. iBellazz ini (2007) recently mea- 
sured the stellar surface density of NGC 2419 to 200 pc. The 
core and half-mass relaxation times of this GC are 9 and 35 
G yr, respectively. 

Mac kev et al.l (120101) recently measured structural and pho- 
tometric properties of MGC1, from which we have learned 
the following. MGC1 resides at approximately 200 kpc from 
M31, and is therefore the most isolated GC known in the Lo- 
cal Group. It has a V-band luminosity of 4 x 10 5 L© and thus 
a stellar mass of w 10 6 M©, a half-mass radius of « 7.5 pc, 
and an indeterminate tidal radius. Mackey et al. have mea- 
sured the stellar surface density for MGC1 out to an impres- 
sive 900 pc. We can estimate the core and half-mass relax- 
ation times of MGC1 by scaling the relaxation times of NGC 
2419 by the 3/2 power of the ratio of their half-mass radii. 
Doing so yields core and half-mass relaxation times of 2 and 
8 Gyr, respectively. 

In Figure [3] we compare the observed stellar surface den- 
sity profiles of NGC 2419 and MGC1 to our model density 
profile for several values of the dark halo-to-stellar mass ra- 
tio, Mdm/M*. We have fixed r s = 250 pc and hence c = 4 for 
simplicity. Such a low value of c is expected for low mass 
halos that formed at high redshift dNavarro et al.l fl997). Data 
are only shown for R p > /?/,. On scales smaller than roughly 
the half-mass radius our assumptions break down (as demon - 
strated by direct TY-body simulations; Baumg ardt et al.l2 002). 

Over the range 10 < R p < 100 pc the data are consistent 
with the predictions for a pure stellar system; models with 
a massive extended dark halo are strongly disfavored. On 
larger scales deviation between the data and models is ap- 
parent, which may be due to tidal stripping or the ongoing 
assembly of the outer stellar halo. 

Our model assumes that a steady state has been achieved 
in the stellar halo via two-body relaxation effects. MGC1 
has a half-mass relaxation time shorter than the age of the 
Universe, and so our technique can be readily applied to this 
system. Indeed, in the absence of a dark halo the effective 
relaxation time decreases with radius (Figure O and so we 
expect a well-developed stellar halo around MGC1. For mod- 
els with massive dark halos (i.e., large Mdm/M*), our results 
strictly apply only to the inner several half-mass radii — at 
larger radii the effective relaxation time becomes longer than 
the age of the Universe. However, even within 1 < R/Rh < 2 
the data strongly favor models with minimal dark matter ha- 
los. Moreover, the data at larger radii are naturally explained 
by assuming that the stellar halo is fully populated by relax- 
ation effects in the absence of an embedded dark halo. 

Conclusions regarding NGC 2419 must be made with 
greater caution than MGC1, since NGC 2419 has a present 
t r that is longer than the age of the Universe. In this case it 
is less clear that our model should apply at all. We are left 
only with the striking agreement between the density profile 
of the halo of this GC and the model prediction that includes 
no dark halo. This strongly suggests, but does not rigorously 
demonstrate, that NGC 2419 is not surrounded by a massive 
dark halo. 

Our results are consistent with Baumgar dt et al.l ((2009), 
who concluded that if a dark halo surrounds NGC 2419, it 
cannot be more massive than 10 7 M© (this is equivalent to a 



limit of Mdm/M* < 10 for this GC). These latter results were 
based on the measured velocity dispersion profile of NGC 
24 1 9 over the range 1 < R p < 60 pc. 

3. DISCUSSION 

In the previous section we argued that the observed stellar 
surface density profiles of the GCs NGC 2419 and MGC1 
place strong constraints on the existence of extended dark 
halos surrounding these GCs. The data are consistent with 
no dark halo, and a firm upper limit on the dark halo mass- 
to-stellar mass ratio is Mdm/M* < 1. The conclusions are 
strongest for MGC1 because it, unlike NGC 2419, has a core 
relaxation time much less than the age of the Universe. 

This upper limit effectively rules out the possibility that 
these GCs formed at the center of their own dark halos, un- 
der the assumption that these GCs have evolved in weak tidal 
fields throughout their lifetimes. This assertion is based on 
the following argument. If these GCs did form within their 
own dark halos and subsequently experienced little tidal strip- 
ping, then the smallest possible value for Mdm/M* would be 
(l-fb)/ft where fb is the universal baryon fraction. Con- 
straint s from the cosmic microwave background imply fb = 
0.17 dKomatsu et al.ll2009l) . and so M DM /M» > 5. Of course, 
less than 100% star formation efficiency, which is expected, 
would only increase this lower limit. Our upper limit of 
Mdm/M* < 1 therefore strongly suggests that these GCs did 
not form within their own dark halos. 

Observations of the outer stellar profile of isolated GCs are 
very sensitive to a dark halo because a dark halo, were it to 
exist, should have a half-mass radius much larger than the 
GC stellar half-mass radius. This fact also explains why it has 
historically been so difficult to obtain strong constraints on the 
presence of a dark halo with kinematic da ta, even with dat a 
extending to several tens of parsecs (e.g.. lLane et al.ll2010h . 
An NFW dark matter halo with a virial mass of 10 8 M© has a 
mass of only 10 6 M© within 50 pc, assuming c = 2 (or within 
10 pc assuming c = 20). For NGC 2419, which has a stellar 
mass of 10 6 M©, the presence of such a halo would be very 
difficult to distinguish from the uncertain corrections required 
to account for low mass stars and stellar remnants, based on 
data that only extends to several tens of pc. 

In recent years it has become clear that most, if not all GCs 
harbor internal spreads in the abundanc e of light elements, 
including CNO, Na, Mg, and Al (see iGratton et ail 12004 
for a review). Several authors have appealed to GC forma- 
tion at the center of extende d dark halos to account for these 
puzzl i ng observations (e.g.. iFreemanl 11993b iBekki & Norrisl 
2006; Bek ki et all 120071; lBokerj|2008t ICarretta et alj|2010alh 
One of the advantages of forming GCs at the center of mas- 
sive dark halos is that they are much less susceptible to ram 
pressure stripping, and, the argument goes, are therefore bet- 
ter able to retain the gaseous material necessary to account 
fo r the observed inte rnal abundance spreads. As discussed 
in IConrov & Sperge| d201 ll) . this line of reasoning is likely 
incorrect because the formation environments of the ancient 
GCs differed substantially from their present day environ- 
ment. The results in this work provide strong independent 
confirmation that indeed GCs which harbor multiple stellar 
populations do not (or need not) form within extended dark 
halos. 

While the current evidence disfavors typical GCs from hav- 
ing formed at the center of their own dark halos, there is some 
reason to suspect that perhaps some of the most massive GCs 
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did indeed form in this way. M54 is the most striking ex- 
ample, as it resides at the center of the disrupting Sagittar- 
ius galaxy, and will in the future likely orbit freely through 
the Galaxy (although recent evidence suggets that M54 re- 
sides at the center of Sagittarius becau se of dynamical fric- 
tion, not because it formed there; see iBellazzini et aDl2008L 
for details). Other candidates for this formation mechanism 
include wCen, M22, NGC 1851, and Gl in M31, all of which 
show internal spreads in the Fe-peak elements. These GCs 
must have formed in deep potential wells in order to retain 
the Fe generated from type la SNe. Nuclear star clusters may 
be the precursors of these massive GCs. The most massive 
GCs in external galaxie s also appear to be self- enriched in Fe 
dStrader & Smith! I2008t iBailin & Harrisl 12009b . although the 
fact that their photometric properties join seamlessly with the 
less massive clusters suggests that GCs of all masses share a 
common origin unrelated to dark halos. Detailed simulations 
will be necessary to conclude whether or not GCs can self- 
e nrich in SNe products w ithout surrounding dark halos. 

Olszewsk i et al.l (120091) recently reported the discovery of 
a 500 pc stellar halo surrounding the GC NGC 1851. Over 
the projected radial range of 50-250 pc, these authors find a 
projected stellar density profile of £ oc r" L24±0 66 . This mea- 
sured profile agrees remarkably well with models that include 
a massive dark halo ( M^m/M* > 10 2 ), which predict a loga- 
rithmic slope of -1.4 over the same radial range. NGC 1851 
currently r esides only 17 kpc from the Galactic center and, ac- 
cording to Olszewsk i et al.ld2009l) . has a period of 0.4 Gyr and 
a perigalacticon of only 5 kpc. The interpretation of the den- 
sity profile of weakly-bound stars in this cluster is therefore 
greatly complicated by the stronger tidal fields it experiences 
and the effect of disk shocking as it crosses the MW disk five 
times per Gyr. The lack of any tidal tails is also peculiar given 
its orbit. As noted above, NGC 1 85 1 shows evidenc e for an 
internal spread in Fe abundance (Carrett a"et alj201 0b). and so 
is a potential candidate for being the remnant of a disrupted 
dwarf galaxy. Future work on the orbit and stellar population 
of this cluster may reveal important clues regarding its forma- 
tion. Radial velocity measurements would be especially valu- 
able, as they should be able to distinguish between a stellar 
halo formed from tidal effects and one formed from loosely 



bound stars on radial orb its. 

Recently, ICohen et al.l ((2010) measured iron and calcium 
abundances of stars in NGC 2419. These authors report the 
discovery of an internal spread in Ca abundances in this clus- 
ter, but no spread in Fe. If confirmed, this result suggests 
that NGC 2419 was able to retain type II SNe ejecta, which 
is difficult to understand unless this cluster was once embed- 
ded within a much deeper potential well than it is currently. It 
could of course be the case that the stars in NGC 2419 simply 
formed from a chemically heterogenous molecular cloud, or 
that the cluster contained many more stars at birth. As with 
NGC 1851, future work on the abundance variations of the 
stars within NGC 2419 and a detailed analysis of its orbit will 
provide essential clues into the origin of this puzzling GC. 

We conclude by recalling a central assumption in the 
present work: that outer halo GCs have evolved in isolation 
throughout their lifetimes. Unless these GCs formed in inter- 
galactic space, they likely once resided within larger proto- 
galactic fragments that have since been tidally destroyed. We 
can say little with confidence regarding the influence of the 
birth environment on the structure of these GCs. Mass lost 
from these young GCs during their first ~ 1 Gyr of evolution 
would result in an expansion of the system due to the loss 
of bin ding energy dKroupa & Boilvl [20021 iMarks & Kroupal 
2010). These effects complicate the interpretation of the outer 
stellar halo of the GCs NGC 2419 and MGC1. Nonetheless, 
the tenuous nature of their stellar halos and the observed sim- 
ilarity in their radial profiles strongly suggests that they are 
being continuously populated by two-body relaxation effects. 
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